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Abstract 

There is a growing interest in producing intense, coherent x-ray radiation with an adjustable 
and arbitrary polarization state. In this paper, we study the crossed undulator scheme (K.-J. 
Kim, Nucl. lustrum. Methods A 445, 329 (2000)) for rapid polarization control in a self-amplified 
spontaneous emission (SASE) free electron laser (FEL). Because a SASE source is a temporally 
chaotic light, we perform a statistical analysis on the state of polarization using FEL theory and 
simulations. We show that by adding a small phase shifter and a short (about 1.3 times the FEL 
power gain length), 90° rotated planar undulator after the main SASE planar undulator, one can 
obtain circularly polarized light - with over 80% polarization - near the FEL saturation. 

PACS numbers: 41.60.Cr 
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I. INTRODUCTION 



Several x-ray free electron lasers (FELs) based on self- amplified spontaneous emission 
(SASE) are being developed worldwide as next-generation light sources . In the soft 

x-ray wavelength region, polarization control (from linear to circular) is highly desirable in 
studying ultrafast magentic phenomena and material science. The x-ray PEL is normally 
linearly polarized based on planar undulators. Variable polarization could in principle be 
provided by employing an APPLE-type undulator {4]. However, its mechanical tolerance 
for lasing at x-ray wavelengths has not been demonstrated, and its focusing property may 
change significantly when its polarization is altered. An alternative approach for polarization 
control is the so-called "crossed undulator" (or "crossed-planar undulator"), which is the 
subject of this paper. 

The crossed-planar undulator was proposed by K.-J Kim to generate arbitrarily polar- 
ized light in synchrotron radiation \s] and PEL sources 6|. It is based on the interference 
of horizontal and vertical radiation fields generated by two adjacent planar undulators in a 
crossed configuration (see Pig. [T]). A phase shifter between the undulators is used to delay 
the electron beam and hence to control the final polarization state. Por incoherent radiation 
sources, the radiation pulses generated in two adjacent undulators by each electron do not 
overlap in time. Thus, a monochromator after the second undulator is required to stretch 
both pulses temporally in order to achieve interference. The degree of polarization is limited 
by beam emittance, energy spread, and the finite resolution of the monochromator, as stud- 
ied in a series of experiments at BESSY 0, ^. On the other hand, for completely coherent 
radiation sources (such as generated from a seeded PEL amplifier or an PEL oscillator), 
the interference occurs due to the overlap of two radiation components in the second un- 
dulator A recent crossed-undulator experiment at the Duke storage ring PEL reported 
controllable polarization switches with a nearly 100% total degree of polarization j^. 

It is well-known in the PEL community that SASE light is transversely coherent but 
temporally chaotic due to the shot noise startup. Thus, the effectiveness of the crossed 
undulator for polarization control deserves a detailed study. In this paper, starting with 
one-dimensional (ID) PEL theory, we calculate both radiation components and generalize 
the results of Ref. [6] to the case of SASE. We then determine the required length of the 
second undulator in order to produce the same average power as that produced in the 
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FIG. 1: (color) Schematic of the crossed undulator for polarization control 

first undulator. We show that the degree of polarization can be determined by the time 
correlation of the two radiation fields and compute its asymptotic expression in the high- 
gain limit. The analytical results are compared with ID SASE simulations after a proper 
statistical averaging. Finally, three-dimensional (3D) effects and simulation results are also 
discussed. 

II. FIELD CALCULATION 

Figure [H shows a schematic of the crossed undulator applied to a SASE FEL. In the 
first planar undulator with a total length Li, spontaneous radiation is amplified to generate 
horizontally polarized SASE field E^. In the second undulator (of length L2) that is rotated 
90° with respect to the first one, propagates freely without interacting with the electron 
beam, while a vertically polarized radiation field Ey is produced by the micro-bunched beam. 
A simple phase shifter such as a four-dipole chicane placing between the two undulators can 
slightly delay the electrons in order to adjust the relative phase of the two polarization 
components. 

In this section, we determine both SASE field components generated by the crossed 
undulator. Let E{z, t) be the complex but slowly varying electric field at undulator distance 
z and time t. We write 

E{z,t) = [ i^^^(^)e*A.[(fc,+fc„).-.ii] ^ 
J V 27r 

where loi = kic is the fundamental resonant frequency corresponding to the average beam 



energy (c is the speed of light); v = uj/uji and Az/ = z/ — 1 is the re^ 



ku = ^tt/Xu with Am the undulator period. Following Refs. 10|, lUl, the ID FEL interaction 



ative frequency detuning. 



starting from shot noise can be described by the coupled Maxwell-Klimontovich equations. 



In the small signal regime before FEL saturation, the equations can be linearized and solved 
by the Laplace transformation: 



where 



J ^TTZ 

F(z) = (f '^^ c-^'P^^"^ KiE,^^dV/d7] - F,(0) 
J 27ri {ri/p - /i) 



2pkuD{jj,) \ 2pku J v/P~l^ 

Di,)=,-^- ji,^f^,. (3) 



2p 7 {r]/p-pf 

Here E'jy and are, respectively, the Fourier components of the electric field and of the 
Klimontovich distribution function that describes the discrete electrons in longitudinal phase 
space, with Ey{Q) and -Fi/(0) the Fourier components of the initial conditions; D{p) = 
determines the FEL dispersion relation where p is the Laplace parameter. In addition, 
parameter p is the dimensionless FEL Pierce parameter V{r]) is the electron energy 
distribution with 7] the relative energy deviation, uq is the electron volume density; ni = 
eK[JJ]/ {A^JqUic^), K2 = e/r[JJ]/(2eo7o), where K is the dimensionless undulator strength 
parameter, the Bessel function factor [JJ] is equal to [Jo(0 ~ "^1(0] with ^ = / {4 + 2K'^) , 
7o is the initial electron energy in units of mc^, and eo is the vacuum permittivity. Note 
that the contour integration of p in Eq. ([2]) must enclose all singularities in the complex p 
plane. Based on this solution, we can calculate radiation field components in the crossed 
undulator according to their initial conditions. 



A. Horizontal radiation field 



The radiation field E^ in the first undulator develops from electron shot noise, with the 
initial conditions 

E:iO) = , / F:iO)dv = TT E ^'^^'''^'^ ' (4) 

where Nx is the number of electrons in one radiation wavelength, and tj{0) is the random 
arrival time of the j*'^ electron at the entrance to the first undulator. We assume the first 



4 



undulator operates in the exponential growth regime. In this regime, the dispersion relation 
has a solution /io with a positive imaginary part that gives rise to an exponentially growing 
field amplitude. For a cold beam with vanishing energy spread, we take V{ri) = 6{ri) in 
Eq. ([2]) and obtain 

In this high-gain regime, the electron distribution from Eq. ([2]) can be simplified as [6|: 

This electron distribution function will be used as an initial condition for the calculation of 
the vertical radiation field as follows. 



B. Vertical radiation field 

The radiation field Ey in the second undulator is generated by the pre-bunched electron 
beam in the first undulator. To control the radiation polarization, the required path length 
delay of the phase shifter chicane is on the order of the PEL wavelength. Such a weak 
chicane does not have dispersive effects that could result in micro-bunching, such as can be 
found for example in an optical klystron (see, e.g., 
at the entrance of the second undulator is 

Ey{o) = o, F,^(o) = . (7) 

As the electron beam develops micro-bunching during the PEL interaction in the first 
undulator, it will radiate coherently in the second undulator. Prom discussions in Ref. 6| 
and simulation results shown in Sec. IIVI below, the intensity of Ey can increase to the same 
level as that of E^ in about one gain length. Thus, for a relatively short second undulator, 
we consider only coherent radiation and ignore any feedback of the radiation on the electron 
beam. With this approximation, the third term at the right hand side of D{fi) in Eq. ([3]) 
can be dropped, and Eq. ([2]) can now be written as 



El{z2) =e 



27ri /i — Az//2p 



arj- 



"^pku J r]/p-iJi 



^ 2mp- Az//2p pI p^ ' ^ ' 
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Here Z2 is the undulator distance from the beginning of the second undulator. The extra 
phase factor e**^ is introduced by the phase shifter just before the second undulator. In the 
last step of Eq. ([HD, we have taken a cold beam with vanishing energy spread and made 
use of the relation niK,2no = 4A;^p^. Note that fiQ is the exponential growth solution that 
satisfies -D(/io) = and is a function of the detuning parameter Au, i.e., 



1^0 ~ -- 



1 - 



Au (Au) 



3p 36p^ 



+ i 



1 - 



{Auf 
36p2 _ 



Eq. (IH]) can be solved by the residue theorem: 



Eyiz2) = E:{L,)e'^*-^'^hmc ( ^ [pKz^ - poe'"(pA;„Z2)'] 



where sinc(x) = sin(x)/a:, ip = Ah'kuZ2, and 



a = arctan 



sin(^/2) 



(9) 



(10) 



(11) 



sinc(-?/'/2) — cos('?/'/2) 

Note that a = 7c/2 when Au = 0. The first term in the square bracket of Eq. (ITOl) describes 
coherent spontaneous radiation from a density-modulated beam and grows linearly with 
the undulator distance Z2 (as discussed in Ref. [l^ in the context of harmonic generation). 
Since the electron beam from the first undulator possesses not only density modulation but 
also energy modulation, the momentum compaction of the second undulator can convert 
the energy modulation into additional density modulation. Thus, the second term in the 
square bracket of Eq. ffTOl) describes the enhanced radiation due to the evolution of the den- 
sity modulations inside the second undulator which grows quadratically with the undulator 
distance. 

In order to generate circularly polarized light, we require that both and Ey have the 
same average amphtude. From Eq. flTUl) . this corresponds to the condition 

2i 



[pkuZ2 - fJ'oe'°'{pkuZ2f] 



1 



(12) 



We consider a cold electron beam with vanishing energy spread, hence the growth rate 
Im(/io) is maximized on resonance, i.e., Au = 0. In this case we obtain the required length 
of the second undulator from Eq. flT^ 



1.3L 



G 



where L 



G 



(13) 



is the ID power gain length. 



III. DEGREE OF POLARIZATION 



E,{t)= j '^Kiz = Li)e^^'^[(^i+'=")^^+'=^^2-'^i*] 



The interference of the two radiation components generated by the crossed undulator will 
produce flexible polarization. At the end of the second undulator when z = Li + L2, these 
radiation fields in the time domain are 

J V 27r 

(14) 

Note that we only used Eq. ([1]) for i^^, at 2; = Li (and ti) and applied the free space 
propagation phase factor e«^!^['=i-t'2-'^i{t-ti)] the second undulator as does not interact 
with the electron beam there. Because of the chaotic nature of SASE radiation, we perform 
a statistical analysis to quantify the state of polarization. 

Following the standard optics textbooks (see, e.g., Refs. jisl, Q), the state of polarization 
can be described by the coherency matrix 



{Eyit)E*it)) {Ey{t)El{t)) 



(15) 



where * means complex conjugate, and the angular bracket refers to the ensemble average. 
The degree of polarization can be calculated as Q 



P= 4 1-4 



detiJl 



(16) 



Xtr[J])2' 

where det[J] and tr[J] are the determinant and trace of the coherency matrix, respectively. 
It is also convenient to introduce the first-order time correlation between E^ and Ey as 

{ES)E;{t)) 



9xy — 



[{mm{\Ey{mYi^ 



(17) 



For polarization control in the crossed undulator, we are particularly interested in the 
case when the average intensities of the two radiation components are the same: {\Ex[t)\^) = 
{\Ey{t)\'^) = I. Under this condition, the coherency matrix simplifies to 



1 

Idxy 1 6 



\9xy\e 



id 



(18) 
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where 9 is the phase difference between and Ey. When ^ = the combined radiation 
is circularly polarized; when ^ = or tt, it is linearly polarized at ±45° relative to the 
horizontal axis. The state of polarization is controllable by adjusting the phase shift (j) in 
Eq. ( ITOl) so that the net phase in g^y is 6* = ±| or O/tt. With equal intensity in both 
transverse directions, the degree of polarization in Eq. (JT6i) is simply given by the amplitude 
of the x-y time correlation, i.e., 

P = \9.y\ ■ (19) 
In the x-ray wavelength region, the electron bunch duration is typically much longer than 
the coherence time of the SASE radiation. Thus, a SASE pulse consists of many random 
intensity spikes that are statistically independent. For a flattop current distribution (of 
width T), we can convert the ensemble average of Eq. (ITTj) into a time average as 

r.T/2 



9xy ^ 



lim i / dtE,{t)E;{t) 
= / u,duE:iL,)Ey*{L,)e-^^-'''-'^\ (20) 

^ J — oo 

where we have applied Eq. f|T^ and the Parseval relation in converting the time integration 
to the frequency integration. Assuming that the first undulator operates in the exponential 
gain regime, the frequency dependence of E^ is approximately Gaussian, i.e., 

IT 

(|i?.^Wr) = ^^e-^, (21) 



where the relative rms SASE bandwidth is 



.27ro:,, 
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III 



Since the short second undulator generates coherent radiation from a pre-bunched beam 
that possesses the same narrow bandwidth (jy, we can expand /ig in Eq. ffTOj) to first order 
in Az/ by using Eq. ([9]). We also ignore the frequency dependence of the second term in the 
square bracket of Eq. flTUl) because its contribution to the radiation intensity is relatively 
small. Finally, we have 

exp (^-^ - i2£^fL2^ sine 



\9xy\ 



oo 

dv- 



(23) 



l + (-2+^-2^3p 

where v = Au/au. In view of Eq. (fT3|) . we take L2 = 1.3i^G in Eq- (|23!) and obtain the 
degree of polarization by computing \gxy\- 
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TABLE I: Main parameters for the LCLS soft x-ray FEL used in simulations. 



Parameter 


value 


unit 


electron beam energy 


4.3 


GeV 


relative energy spread 


0(0.023) 


% 


bunch peak current 


2 


kA 


transverse norm, emittance 


1.2 


/im 


average beta function 


8 


m 


undulator period 


3 


cm 


undulator parameter K 


3.5 




FEL wavelength 


1.509 


nm 


FEL p parameter 


0.119 


% 


ID power gain length Lq 


1.17 


m 


3D power gain length Lq^ 


1.48 


m 



IV. NUMERICAL SIMULATIONS 



A. ID results 



We first use a ID FEL code to simulate the SASE radiation produced by the crossed 
undulator configuration and to analyze the degree of polarization. The code follows the 
time-dependent approach developed in Ref. 17| and employs the shot noise algorithm of 
Penman and McNeil Electron energy spread can be included using Fawley's beamlet 



method 



19| . After computing the field produced in the first undulator, we allow to 



propagate freely without further interacting with the electron beam. The simulated electron 
distribution from the first undulator is then used to generate the Ey field in the second 
undulator. 

As a numerical example, we use the parameter set listed in Table [T] that is similar to the 
soft x-ray LCLS operation In the ID simulations, the energy spread is set to zero since 
we want to compare with the previous analytical results. 

Fig. [2] shows the average radiation power in both x and y directions produced by the 
cross undulator. The length of the first undulator is allowed to vary, while the second 
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FIG. 2: (color) ID simulations of the average SASE power at 1.5 nm from the first (blue cross) and 
the second (red plus) undulator. Here Li is the length of the first undulator, L2 = l.SLc = 1.53 
m is the length of the second undulator. 

undulator length L2 = I-SLq ~ 1.53 m is held constant. As predicted by Eq. ( |T3l) . the 
power of the two radiation components are essentially the same in the exponential gain 
regime. Near saturation, the power of the vertical field is lower than that of the horizontal 
one because the FEL-induced energy spread starts to de-bunch the electron beam in the 
second undulator. We repeat the simulations 200 times for each Li with different random 
seeds to start the process and calculate the first-order time correlation between and Ey 
at the exit of the second undulator using the ensemble average defined in Eq. ( IT71) . Figure [3] 
shows the amplitude of this correlation from the simulation results as well as the numerical 
integration of Eq. fl23l) (the red solid curve) for a comparison. When the first undulator 
is less than a couple of gain lengths, the crossed undulator operates in the spontaneous 
emission regime, the amplitude of the x-y correlation and hence the degree of polarization 
are very small without the use of a monochromator. The degree of polarization increases in 
the exponential growth regime and reaches a maximum of 85% near the FEL saturation. In 
this regime and especially when the gain is very high, we see very good agreement between 
simulations and Eq. (123|) . In the saturation regime, the amplitude of the x-y correlation 
starts to decrease, and the linear theory starts to deviate from the simulation results. 

There are two effects that prevent the degree of polarization to reach 100% in a crossed- 
undulator SASE FEL. First, there is relative slippage between E^ and Ey in the second 
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FIG. 3: (color) Amplitude of the time correlation \gxy\ from the ID simulations (black star) and 
from Eq. ([23]) (red solid curve). The degree of polarization is equal to \gxy\ when (P^.) = {Px)- For 
comparison, the blue dashed curve shows an estimate of the polarization by Eq. (j24p (see text for 
more details). 

undulator. Since Ex stops interacting with the electron beam after the first undulator, the 
group velocity of is the speed of light c. However, the group velocity of Ey is slower than 
c because it is generated by the micro-bunched beam that travels at the average longitudinal 
velocity (3\\c. In fact, ID simulations indicate that the group velocity of Ey is almost the 
same as that of the electrons within the short second undulator section. (This numerical 
result is also confirmed in 3D simulations to be discussed in the next section.) To estimate 
the slippage effect, we take Ey{t) ^ E^it — r) with cr = ^2(1 — /3||) and apply the first-order 
time correlation function of the SASE field to estimate \gxy\'- 

-2 • 



g{T) = exp 



vrr 
'2^ 



(24) 



where Tc = yvr/c"w is the coherence time [20[, and a^^ is given by Eq. ( l22l) . Equation (E 



yields the blue dashed curve shown in Fig. [3l which indicates that the slippage effect only 
accounts for about a half of the depolarization in the crossed undulator. A careful exami- 
nation of the intensity profile between Ex and Ey shows a visible difference from a simple 
time delay (see Fig. Hlfor a 3D example). This accounts for the additional depolarization 
effect in a crossed undulator SASE FEL. 
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B. 3D Discussions 



A remarkable feature of a SASE FEL is its transverse coherence. At a sufficiently high 
gain, a single transverse mode with the largest growth rate will dominate over all other 
transverse modes for a typical SASE FEL. Thus, we expect the previous ID analysis still 
applies to 3D situations in the high gain limit, with the maximum polarization obtainable at 
the end of the exponential growth regime. Since the length of the second undulator is short, 
the diffraction effects for the free-propagating in the x-ray wavelength regime is expected 
to be small. Thus, the 3D effects such as emittance and diffraction do not play significant 



roles in determining the degree of polarization 



We use the 3D FEL code GENESIS L3 



br a crossed undulator SASE FEL. 



2l| to check these expectations. The electron 



beam is dumped at the end of the first undulator and is used to generate Ey in the second 
undulator. E^. propagates in the same length of the second undulator but without any 
undulator magnetic field. We use the same soft x-ray FEL example listed in Table [T] as the 
ID case but with a relative energy spread of 0.023%, which roughly corresponds to the LCLS 
soft x-ray parameters. The length of the first undulator is chosen to be 23 m and is about 3 m 
before the saturation point. A 2-m short second undulator is necessary to produce the same 
radiation power for the vertical field (see Fig.H]). The 3D power gain length corresponding to 
these parameters is L^P = 1.48 m, so Eq. f|T3|) approximately holds in this 3D case. We use 
the on-axis far-field radiation intensity and phase from GENESIS simulations to calculate 
the time correlation between Ex and Ey of Eq. ( |T7I) . Instead of performing many statistical 
runs for the ensemble average, we average the result over hundreds of intensity spikes within 
the radiation pulse in order to save on simulation effort. The amplitude of the x-y correlation 
from this 3D calculation is 87%, very close to the ID prediction. Figure H] shows the central 
section of the simulated power profiles and Py at the end of the second undulator. A 
small time delay due to the slippage effect and a somewhat different temporal structures 
between P^ and Py are the main depolarization effects, as discussed in the previous section. 

V. CONCLUSIONS 

The statistical analysis presented in this paper shows that the crossed-planar undulator 
is an effective method for polarization control in a SASE FEL. To optimize the degree of 
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FIG. 4: (color) GENESIS simulated power profiles of the horizontal field (blue dashed curve) and 
the vertical field (red solid curve) at the end of the second undulator, bunch head at left. 

polarization, the first undulator should operate at the end of the exponential growth regime, 
while in order to generate circularly polarized x-rays, the second undulator should be about 
1.3 times the power gain length. The maximum degree of polarization is over 80% from both 
theory and simulations. If fast pulsed magnets are employed in the phase shifter chicane, the 
relative phase between the two radiation components from the crossed undulator can vary at 
hundreds of Hz, hence enabling fast polarization switching for many scientific applications. 
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